
HUAPIU_TS = reshape(HUAPIU',[1,1739*1440]);

Hd=myfilter;

HUAPIU_FILT_TS = filter(Hd,HUAPIU_TS);

HUAPIU_FILT_TS(1:653+1440) = []; %first data is 2000/8/26/ 00:00:30
HUAPIU_FILT_TS(end-786:end) = [];%last day is 2005/12/12 23:59:30
HUAPIU_FILT = reshape(HUAPIU_FILT_TS,[1440,1737])';
HUAPIU_FILT_fday = HUA_PIU_fday(2:end-1);


%fday_huapiu_filt = (datenum(2000,8,27,0,0,30):(60/(24*3600)):datenum(2005,12,11,23,59,30))-datenum(2000,1,1);




ESKH = reshape(ESK(:,5:end)',[24*60,2192]);
L = ESKH>99999;
ESKH(L) = NaN;
fday_esk = [0:2191];
base = repmat(mean(ESKH(60:160,:))',[1,1440])';
ESKH_dH = ESKH-base;



